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Abstract — Homology theory provides new and powerful so- 
lutions to address the coverage problems in wireless sensor 
networks (WSNs). They are based on algebraic objects, such 
as Cech complex and Rips complex. Cech complex gives accu- 
rate information about coverage quality but requires a precise 
knowledge of the relative locations of nodes. This assumption is 
rather strong and hard to implement in practical deployments. 
Rips complex provides an approximation of Cech complex. It is 
easier to build and does not require any knowledge of nodes 
location. This simplicity is at the expense of accuracy. Rips 
complex can not always detect all coverage holes. It is then 
necessary to evaluate its accuracy. This work proposes to use 
the proportion of the area of undiscovered coverage holes as 
performance criteria. Investigations show that it depends on 
the ratio between communication and sensing radii of a sensor. 
Closed-form expressions for lower and upper bounds of the 
accuracy are also derived. For those coverage holes which can 
be discovered by Rips complex, a homology-based distributed 
algorithm is proposed to detect them. Simulation results are 
consistent with the proposed analytical lower bound, with a 
maximum difference of 0.5%. Upper bound performance depends 
on the ratio of communication and sensing radii. Simulations also 
show that the algorithm can localize about 99% coverage holes 
in about 99% cases. 

Index Terms — Wireless sensor networks, coverage hole, homol- 
ogy- 

I. Introduction 

WIRELESS sensor networks (WSNs) have attracted a 
great deal of research attention due to their wide poten- 
tial applications such as battlefield surveillance, environmental 
monitoring and intrusion detection. Many of these applications 
require a reliable detection of specified events. Such require- 
ment can be guaranteed only if the target field monitored by 
a WSN contains no coverage holes, that is to say regions of 
the domain not monitored by any sensor Coverage holes can 
be formed for many reasons, such as random deployment, 
energy depletion or destruction of sensors. Consequently, it 
is essential to detect and localize coverage holes in order to 
ensure the full operability of a WSN. 

There is already an extensive literature about the cover- 
age problems in WSNs. Several approaches are based on 
computational geometry with tools such as Voronoi diagram 
and Delaunay triangulations, to discover coverage holes ([T|- 
These methods require precise information about sensor 
locations. This substantially limits their applicability since 
acquiring accurate location information is either expensive or 
impractical in many settings. Some other approaches attempt 
to discover coverage holes by using only relative distances 
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between neighbouring sensors p^-fSI. Similarly, obtaining 
precise range between neighbour nodes is costly. 

More recently, homology is utilized in |j7)-|[9) to address the 
coverage problems in WSNs. Christ and his collaborators in- 
troduced a combinatorial object, Cech complex (also known as 
nerve complex), which fully characterizes coverage properties 
of a WSN (existence and locations of holes). Unfortunately, 
this object is very difficult to construct as it requires rather 
precise information about the relative locations of sensors. 
Thus, they introduced a more easily computable complex. Rips 
complex (also known as Vietoris-Rips complex). This complex 
is constructed with the sole knowledge of the connectivity 
graph of the network and gives an approximate coverage by 
simple algebraic calculations. As regards implementation in 
real WSN, these homology based methods are necessarily 
centralized, which makes them impractical in large scale 
sensor networks. Some algorithms have been proposed to 
implement the above mentioned ideas in a distributed context, 
TOl, |11|. But there are two disadvantages of these 
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algorithms. On one hand, these homology based algorithms 
are all dependent on the assumption that the communication 
radius of a sensor is smaller than \/3 times the sensing radius 
of the sensor When such assumption is not satisfied, it is 
possible that Rips complex may miss some special coverage 
holes (such holes are defined as triangular holes in Section 
III). It is thus of paramount importance to determine the 
proportion of missed coverage holes to assess the accuracy 
of Rips complex based coverage hole detection. On the other 
hand, these algorithms are either slow or require coverage 
holes separated with enough distance. So it is necessary to 
design an efficient and accurate algorithm to detect coverage 
holes. 

The main contributions of our paper are as follows. First, 
The relationship between Cech complex and Rips complex 
in terms of coverage hole is analysed under different ratios 
between communication and sensing radii of a sensor We 
find that when the communication radius is at least two times 
sensing radius, if there is a hole in Rips complex, there must 
be a hole in Cech complex. A hole in a Cech complex missed 
by a Rips complex must be bounded by a triangle. Based on 
that, a formal definition of triangular and non-triangular hole 
is presented. 

Second, the proportion of the area of triangular holes under 
a homogeneous setting is analysed. It is indicated that such 
proportion is related to the ratio between communication 
and sensing radii of each sensor In addition, closed-form 
expressions for lower and upper bounds of the proportion 
under different ratios between communication and sensing 
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radii are derived. 

Third, a homology based distributed algorithm is proposed 
to detect non-triangular coverage holes. In the algorithm, a 
Rips complex is first constructed for a given WSN. Then some 
vertices and edges are deleted without changing the number 
of holes in the original Rips complex. After that, the edges 
lying on the boundary of holes will be detected. Then coarse 
boundary cycles can be discovered. Finally all boundaries of 
the non-triangular holes are found by minimizing the length 
of coarse boundary cycles. 

The remainder of the paper is organized as follows. Section 
II presents the related work. In Section III, the network model 
and the formal definition of triangular and non-triangular hole 
are given. Upper and lower bounds on the proportion of 
the area of triangular holes under different ratios between 
communication and sensing radii are computed in Section IV. 
Section V describes the homology based distributed algorithm 
for non-triangular holes detection. In Section VI, performance 
evaluation of the bounds and the algorithm is given. Finally, 
Section VII concludes the paper 

II. Related Work 

There is not much work on the proportion of the area of 
triangular holes. Some recent work p2)-p3| provides some 
results on coverage probability but with a different point of 
view. In | [T2[ , | [T3|, t he fraction of the area covered by sensors 
is analysed. In Jl4[, the authors study how the probability 
of coverage changes with the sensing radius or the number 
of sensors. In p3] , a point in a plane is defined to be tri- 
covered if it lies inside a triangle formed by three nodes, 
and the probability of tri-coverage is analysed. None of them 
considers triangular holes, we provide some initial results 
about proportion of triangular holes in |jT6j. 

The related work on the coverage problem in WSNs can 
be generally classified into three categories: location-based, 
range-based and connectivity-based. In the first category, the 
Voronoi diagram is used to detect boundary nodes in |[T|, Q. 
In Q, a boundary node detection algorithm based on localized 
Voronoi polygons is proposed. In p7) , it is proved that a 
sensor node does not border a coverage hole if its sensing 
border is entirely covered by the sensing ranges of its neigh- 
bours. Another boundary node detection approach proposed 
in p8) simplifies the previous border-checking approach by 
only checking intersection points on the sensing border The 
former algorithms only detect boundary nodes but do not 
detect boundary cycles. Based on that, the authors propose a 
method to obtain the exact boundary cycles of coverage holes 
in |19|. 

Range-based approaches attempt to identify boundary nodes 
based on relative distance between neighbouring nodes. In 
Q, a localized Voronoi polygon based boundary node de- 
tection algorithm is proposed, which is similar as that in 
p. The difference lies in that the localized Voronoi polygon 
is constructed using location information of nodes in |j3| 
while it is constructed based on directional and distance 
information between neighbouring nodes in Q. In |j5), ||6), 
the author proposes a coverage verification algorithm based on 



distances between neighbouring nodes. These approaches also 
only detect boundary nodes without discovering all boundary 
cycles. 

Connectivity-based schemes try to discover boundary cycles 
using only connectivity information. Our algorithm belongs to 
this category. As a pioneer work, in I^TJ, Ghrist et al propose 
an algorithm that detects coverage holes via homology. They 
construct the Rips complex corresponding to the commu- 
nication graph of the network and determine the coverage 
by verifying whether the first homology group of the Rips 
complex is trivial. But their method is centralized. The first 
steps of implementing the above ideas in a distributed way 
are taken in |10|. It is shown that combinatorial Laplacians 
are the right tools for distributed computation of homology 
groups and can be used for decentralized coverage verification. 



In 1 20 1, a gossip-Uke decentralized algorithm for computation 



of homology groups is proposed. But its convergence is slow. 
In 1 11 1, a decentralized scheme based on Laplacian flows is 
proposed to compute a generator of the first homology. And 
they formulate the problem of localizing coverage holes as an 
optimization problem for computing a sparse generator of the 
first homology. But it is possible that some cycle found by 
their algorithm encircles more than one coverage holes. More 
recently, some distributed algorithms are proposed to detect 
topological holes pT[ , [22) and some other algorithms are 
proposed for boundary recognition in sensor networks by using 
only connectivity information p3)-|[26). These algorithms 
may be useful for coverage hole detection, but they can only 
detect large coverage holes and may miss small area coverage 
holes. The algorithm proposed in [27J can detect all non- 
triangular holes, but our algorithm is more efficient, especially 
when coverage holes are largely separated. 

III. Models and Definitions 

A brief introduction to the tools used in the paper is given. 
For further readings, see ||28|-p0|. Given a set of points V, a 
^-simplex is an unordered set [vq,vi, ...,Vk] Q V where Vi ^ 
Vj for all i 7^ j. So a 0-simplex is a vertex, a 1 -simplex is an 
edge and a 2-simplex is a triangle with its interior included, see 
Figure [T] In the following part, a triangle refers to a 2-simplex 
if not indicated explicitly. The faces of this ^-simplex consist 
of all (fc-l)-simplex of the form [wq, Wi-i, Wi+i, Wfc] for 
< i < k. A simplicial complex is a collection of simplices 
which is closed with respect to inclusion of faces. 
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Cech complex and Rips complex are two simplicial com- 
plexes defined as follows Q. 

Definition 1 (Cech complex). Given a collection of sets U, 
Cech complex ofU, C{U), is the abstract simplicial complex 
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whose k-simplices correspond to non-empty intersections of k 
+ 1 distinct elements ofU. 

Definition 2 (Rips complex). Given a set of points X in K" 
and a fixed radius e, the Rips complex of X, TZ^(X), is the 
abstract simplicial complex whose k-simplices correspond to 
unordered (k + l)-tuples of points in X which are pairwise 
within Euclidean distance e of each other 

Consider a collection of stationary sensors (also called 
nodes) deployed randomly in a planar target field. As usual, 
isotropic radio propagation is assumed. Each sensor monitors 
a region within a circle of radius Rs and may communicate 
with other sensors within a circle of radius Re- 
in addition, some other assumptions are as follows. 

1) There are sensors located on the external boundary of 
the target field. They are known as fence sensors and 
other sensors are referred to as internal sensors. Each 
fence sensor has two fence neighbours. This is also the 
general assumption in many homology based algorithms 

il7i-f9i, nn. 

2) Although sensors are not aware of their locations, every 
sensor can know whether it is a fence or an internal 
node by using the mechanisms presented in fSl or other 
methods as in [25 J . In fact, it is a conventional assump- 
tion adopted by many existing range-based methods ||5], 
(|3T1 or connectivity methods 1 11 1, |25|. 

3) Internal sensors are distributed in the planar target field 
according to a homogeneous Poisson point process with 
intensity A. 

4) Each sensor has an unique ID. 

5) The network has only one connected component. 

Let V denote the set of sensor locations in a WSN and 
S = {stj, w e V} denote the collection of sensing ranges of 
these sensors. For a location v, s^, = {x £ : ||a;— w|| < Rs]- 
Then, according to the definition, the Cech complex and 
Rips complex of the WSN, respectively denoted by Cfl^(V) 
and TZn^iy), can be constructed as follows: a fc-simplex 
[vq.vi,--- ,Vk\ belongs to Cr^{V) whenever H^'^^qS^, 7^ 
and a A;-simplex [vq, ui, • • • , v^] belongs to TZj^^iV) whenever 
\\vi - Wm[[< Rc for all < Z < m < fc. 

Figure |2| shows a WSN, its Cech complex and two Rips 
complexes for two different values of Depending on the 
ratio Rc over Rs, the Rips complex and the Cech complex may 
be close or rather different. In this example, for Rc = 2Rs, 
the Rips complex sees the hole surrounded by 2,3,5,6 as in 
the Cech complex whereas it is missed in the Rips complex 
for Rc = 2.5Rs- At the same time, the true coverage hole 
surrounded by 1, 2,6 is missed in both Rips complexes. 

In fact, as proved in any coverage hole can be found 
in Cech complex. Furthermore, there are following relations 
between Cech complex and Rips complex: 

7^i^.(V) c Cr^V) c 7^27^,(V), ifi?, < ViRs- (D 

According to ([I}, some relationships between Cech complex 
and Rips complex in terms of coverage hole can be derived as 
illustrated in the following corollaries. For convenience, define 

7 = Rc/Rs. 




(c) 



(d) 



Fig. 2. (a) a WSN, (b) Cech complex, (c) Rips Complex under Rc = 2Rs, 
(d) Rips Complex under Rc = 2.5-Rs 



Corollary 1. When 7 < if there is no hole in Rips complex 
TZr^{V), there must be no hole in Cech complex CR^iV). 

Proof: If there is no hole in TZR^iV), it means that 
Ti-R^iy) can be triangulated. Since 7 < \/3 means Rc < 
V^Rs, according to the first inclusion in ([T]), we can see that 
TZr^{V) C Cr^ (V). Consequently, Cech complex Cr^ (V) can 
also be triangulated. And when Rc < V^Rs, each triangle 
must be covered by the sensing range of its vertex nodes |7|. 
So there is no hole in Cfl^(V). ■ 

Corollary 2. When "f > 2, if there is a hole in Rips complex 
TZr^ (V), there must be a hole in Cech complex Cr^ (V). 

Proof: If there is a hole in TIr^{V), there must be a 
cycle with more than three edges in TIr^{V) that can not be 
triangulated, as the cycle {2, 3, 5, 6} in Figure ^c)- Since 
7 > 2 means Rc > 2Rs, according to the second inclusion 
in ([1}, we can see that Cr^V) C 7^2i^,(V) C 7^i^,,(V). 
Consequently, there must also be a cycle in Cr^ (V) which 
can not be triangulated. And there is a coverage hole in the 
cycle. ■ 

Corollary 3. When \/3 < 7 < 2, there is no guarantee 
relation between Rips complex TZr^ (V) and Cech complex 
Cr,. (V) in terms of holes- 

Proof: It is a direct corollary from Corollary [T] and |2] ■ 
From the discussion above, a hole in a Cech complex not 
seen in a Rips complex must be bounded by a triangle. Based 
on this observation, a formal definition of 'triangular hole' and 
'non- triangular hole' is given as follows. 

Definition 3 (Triangular and non-triangular hole). For a pair 
of complexes Cr^ (V) and TZr^ (V), a triangular hole is 
an uncovered region bounded by a triangle which appears 
in TZr^{V) but not in Ci^^(V). Any other holes are non- 
triangular 

For triangular holes, it is impossible to detect them with 
only connectivity information, so we want to analyse the 
proportion of the area of such holes in a target field. For non- 
triangular holes, we aim to design a distributed algorithm to 
discover the boundaries of these holes. 
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IV. Bounds on proportion of the area of 

TRIANGULAR HOLES 

In this section, the conditions under which any point in the 
target field is inside a triangular hole are first given. From 
the discussion in Section III, it is found that the proportion 
of the area of triangular holes is related to the ratio 7. Three 
different cases are considered for the proportion computation. 
For each case, the upper and lower bounds of the proportion 
are derived. 

A. Preliminary 

Lemma 1. For any point in the target field, it is inside a 
triangular hole if and only if the following two conditions are 
satisfied: 

1) the distance between the point and its closest node is 
larger than R^. 

2) the point is inside a triangle: the convex hull of three 
nodes, two by two less than or equal to Rf. apart. 

Lemma 2. If there exists a point O which is inside a triangular 
hole, then Rg < Rc/V^- 

Proof According to the definition of triangular holes, if 
there exists a triangular hole, then there must be a triangle A 
satisfying A e 7^J^^(V) and A ^ Cr^{V). If Rs > Rc/VS, 
then according to the first inclusion in (1), we have TZn^iV) G 
Cr^{V), it means that there exists not a triangle which is in 
TZr^{V) but not in C^, (V), there is a contradiction, so Rg < 

rJVs. " ■ 

Lemma 3. Let O be a point inside a triangular hole and I 
denote the distance between O and its closest neighbour, then 
Rs<l< Rc/VS- 

Proof Rs < I is a direct corollary from Lemma [T| We 
only need to prove I < i?c/\/3- If point O is inside a triangular 
hole, it must be surrounded by a triangle formed by sensors 
with pairwise distance less than or equal to Re- Assume it is 
surrounded by a triangle N0N1N2, as in Figure |3] The closest 
neighbour of O is not necessarily in the set {No,Ni,N2}. If 
I > Rc/VS, then dg > I > Rc/V^, di > I > Rc/V^ and 
d2>l> Rc/V^- In addition, since ZNqONi + ZN1ON2 + 
ZNqON2 — 2tt, there must be one angle no smaller than 
2tt/3. Assume ZNqON2 > 27r/3 and denote it as a. Then 
according to the law of cosines, = (iQ+c?| — 2c?od2 cos a > 
R^J3+RI/3-2/3RcRcCos{2tt/3) ^ R^. So do2 > Rc- Since 
No and N2 are neighbours, do2 < Rc- There is a contradiction. 
Therefore I < Rc/V3. ■ 

A Poisson point process whose intensity is proportional 
to the Lebesgue measure is stationary in the sense that any 
translation of its atoms by a fixed vector does not change its 
law. Thus without considering border effect p2) , any point 
has the same probability to be inside a triangular hole as the 
origin O. This probability in a homogeneous setting is also 
equal to the proportion of the area of triangular holes. We 
borrow part of the line of proof from |15| where a similar 
problem is analysed. 

We consider the probability that the origin O is inside a 
triangular hole. Since the length of each edge in the Rips 



Ni 




Fig. 3. Illustration of O being inside a triangular hole 

complex must be at most Rc, only the nodes within Rc 
from the origin can contribute to the triangle which bounds 
a triangular hole containing the origin. Therefore, we only 
need to consider the Poisson process constrained in the closed 
ball B{0,Rc) which is also a homogeneous Poisson process 
with intensity A. We denote this process as $. In addition, 
T{x,y,z) denotes the property that three points x,y,z are 
within pairwise Euclidean distance Rc from each other, and the 
origin O is inside the triangular hole bounded by the triangle 
with these points as vertices. When no, ni, ^2 are points of the 
process T{no,ni,n2) is also used to denote the event that 
the triangle formed by the nodes no, ni, n2 bounds a triangular 
hole containing the origin. In addition, we use T'{nQ,ni,n2) 
to denote the event that the nodes no,ni,n2 can not form a 
triangle which bounds a triangular hole containing the origin. 

Let tq = To{^) be the node in the process $ which is 
closest to the origin. There are two cases for the origin to 
be inside a triangular hole. The first case is that the node tq 
can contribute to a triangle which bounds a triangular hole 
containing the origin. The second case is that the node tq can 
not contribute to any triangle which bounds a triangular hole 
containing the origin but other three nodes can form a triangle 
which bounds a triangular hole containing the origin. So the 
probability that the origin is inside a triangular hole can be 
defined as 

p(A,7) = P{0 is inside a triangular hole} 

= P{ IJ T(no,ni,n2)} 

{no,rii,n2}C<I> 

= P{ U r(ro,ni,n2)} 

{ni,ri2}C$\{To(*)} 

+ P{ IJ T(no,ni,n2) I r'(ro,n3,n4)} 

{no,--- ,n4}C*\{ro(*)} 

(2) 

In the following parts, we will analyse this probability in 
three different cases. 

B. Case < 7 < ^3 

Theorem 1. When < 7 < \/3, p(A, 7) = 0. 

Proof: It is a direct corollary from Lemma |2] ■ 
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C. Case \/3 < 7 < 2 

Theorem 2. When < 7 < 2, pi{X,^) < p(A,7) < 
p„(A,7), where 

rRo/V3 rai /.fli(ro,ei) 

P,(A,7) = 2^A2 / rodro / d9, / 

anc/ 

/•Rc/\/3 rai /•Ri{ro,0i) 

+ P{ IJ r(no, 711,^2) I T'(To,n3, 

{no,-,«4}C$\{ro(*)} 

anc/ 

ao = 2 arccos(i?c/(2ro)) 
tti = 2arcsin(i?c/(2ro)) — 2 arccos(i?c/(2ro)) 



formed by tq, ni and n2, ni must lie in the shadow area 
A+ = H+f]BiTo,Rc)f]B{0,Rc)\BiO,do)f]B{M2,Rc), 
shown in Figure |4] A/2 is one intersection point between the 
circle C{0,do) and the circle C{to,Rc), such intersection 
point must exist in this case since Rc = jRs < 2i?s < 2do. 



H* 







/ \ 


/ / \ 1 A 

^// \ -/ 






^ ^ ^ 'ro 


mrfo, It) ^ ^ ao 1 ; 






\ R< , 
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R\{rQ,d\) = min(A/i?2 — rQ sin^ — rocos^i 



Fig. 4. Illustration of areas A+ , 5+ and S in the case V3 < 7 < 2 



\J Rl — r"^ sin^(0i + ao) + ro cos(0i + ao)) Ordering the nodes in A+ by increasing polar an- 
gle so that Ti = {di^Oi) has the smallest angle 6i. 
And assume the nodes ro, ri and another node T2 G 
i7"n-S(To,i?c)n-B(0,i?c)\B(0,do) can form a triangle 
which bounds a triangular hole containing the origin, then T2 
must lie to the right of the line passing through ti and O, 
denoted by H^(di) which contains all points with polar angle 
9 £ {9i — TT,9i). In addition, the distance to ti is less than 
Rc- So the node T2 must he in the region S^, as illustrated 
in Figure |4] 



Proof: We first prove the lower bound. It can be obtained 
from (|2| that 

p(A,7)>P{ y T(ro,ni,n2)} 

{ni,n2}C*\{To(*)} 



So for the lower bound, we only consider the case that the 
closest node tq must contribute to a triangle which bounds a 
triangular hole containing the origin. 

Using polar coordinates, we assume the closest node tq lies 
on {do,TT). It is well known that the distance do is a random 
variable with distribution 



Fa„ (ro) = P{do < ro} - 1 - g-^^'-o (5) 
Therefore the above probability can be written as 



P{ U T{ro,ni,n2)} 

{ni,n2}C*\{To(*)} 

P{ y r((ro,7r),ni,n2)}Fd„(dro) 

{ni,n2}C*' 



(6) 



where ^J,^^ is the restriction of in B{0,Rc)\B{0,ro). 

Once the node To is determined, the other two nodes must 
lie in the different half spaces: one in = M+ x (0, tt) and 
the other in — M+ x {—it, 0). Assume rii lies in iJ+ and 
n2 lies in . Since the distance to tq is at most R,., rii 
and n2 must also lie in the ball B{tq,Rc). Furthermore, the 
distance to the origin is at most Rc and larger than do, they 
should also lie in the area B{0, Rc)\B{0,do). Therefore, 
ni must lie in H+ f]B{To,Rc)r\BiO,Rc)\B{0,do) and 
n2 must lie in H- f] B{to, R^) f] B{0, Rc)\B{0, do). In 
addition, considering the distance between ni and n2 should 
be at most Rc and the origin should be inside the triangle 



5-(to,ti) = S-ido,d,,9,) = H-f]BiTo,Rc) 
fl 5(0, Rc)\BiO, do) fl H+{9i) fl B{n,Rc) 

Here we need to obtain the density of node ti . Considering 
the way ti was defined, there should be no nodes with a polar 
angle less than 9i, that is to say no nodes are in the region 

S+{To,n) = S+{do,9,) ^ A+ f]H+{9i) 

Since the intensity measure of the Poisson process in polar 
coordinates is \rdrd9, the density F^.-^ of ti can be given as 

Fr,{dri,d9i) = Arie-^l'5^(*^'^i)ldrid6'i (7) 

The integration domain D{do) with respect to parameters 
(dij^i) can be easily obtained. From the construction of the 
area , we can get ao — 2 arccos(_Rc/(2do)) ™d ai = 
2 arcsin(i?c/(2(io)) — 2 arccos(_Rc/(2do))- So ao < 6*1 < ai 
and do < di < Ri{do,9i), where 



Ri{dQ,9i) =mm{\/ — dg sin 9i — docos^i, 



\J Rl-dl sin^(6li + ao) + do co&{9i + ao)) 
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Assume only to,ti and nodes in 5^(to,ti) can contribute 
to the triangle which bounds a triangular hole containing the 
origin, we can get a lower bound of the probability that the 
origin is inside a triangular hole. It is a lower bound because 
it is possible that ti can not contribute to a triangle which 
bounds a triangular hole containing the origin, but some other 
nodes with higher polar angles in the area can contribute 
to such a triangle. For example, in Figure [5] if there is no 
node in but there are some nodes in S'^, then ti can 
not contribute to any triangle which bounds a triangular hole 
containing the origin, but t{ can form such a triangle with tq 
and another node in S'^. Based on the assumption, we have 

P{ U T((ro,^),ni,n2)} 

>P{ U T((ro,7r),ri,n2)} 

P{ U T{{ro,7r),{n,9i),n2)}Fr,idn,d0,) 

S-(ro,ri,ei) 

J JD(ra) 




Fig. 5. Illustration of areas 5"+ and S' in the case < 7 < 2 

Therefore, from Q, (|7]i and (|8]l, the lower bound shown 
in (|3]l can be derived. 

Next we will prove the upper bound. As discussed in 
Section IV.A, there are two cases for the origin being inside 
a triangular hole. As for the second case that the closest node 
To can not but some other nodes can contribute to a triangle 
which bounds a triangular hole containing the origin, it is not 
easy to obtain a closed-form expression for such probability, 
but we can get it by simulations. Simulation results show that 
this probabiUty is less than 0.15% at any 7 < 3 with any 
intensity A. So we still focus on the probability of the first 
case. 

For the lower bound, we only considered the case that 
Ti contributes to a triangle which bounds a triangular hole 
containing the origin. For the upper bound, we need to further 



consider the case that ti can not but some other nodes in 
can contribute to such a triangle, shown in Figure |5] Assume 
the node t[ = {d\,6'^) with the second smallest polar angle in 
can contribute to such a triangle, it means that there is no 
node in {do ,di,6i) but there is at lease one node in the re- 
gion S''ido,d,,di,d[,0[) ^S-{do,d[,0[)\S-{do,di,9i). 
Then the density of the pair {ti,t[) is given as 



Fr,,riidri,dei,dr[,de[) 

= X\ir[e~^\^^'-^°-'''^^\dnd9idr[de[ 



(9) 



The probability that ti can not but t[ can form a triangle 
which bounds a triangular hole containing the origin with tq 
and another node in {do, di, 9i, d'l, 9'i) can be given as 

P{ U T((ro,7r),T{,n2) |r'((ro,7r),ri,n3)} 



P{$;j5-(ro,ri,0i))-O}x 
P{$;„ {S'- (ro ,n,9,,r[,9[))> 0}F,^ {dr, , d9, , dr[ , d9[ ) 

,-\\S-(ro.ri,ei)\^ 

(8) (^i_^-MS'-{ro,rueur[,e[)\s^p^^ ^,{dn, d9i, dr[, d9[) 

(10) 

As we can see from Figure |5] as long as t( has a 
higher polar angle than ri has, the sum of \S^ {ro,r[,9'i)\ 
and \S'^ {ro,ri,9i,r'i,9'i)\ will be always smaller than 

\S-{ro,ro,9,)\. 



Therefore we can get from (lOi 



P{ U T{{rQ,TT),T[,n2) \T'{{ro,TT),Ti,n3)} 



(11) 



Fr,y{dri,d9i,drld9[) 
Furthermore, let S'+{da,9i,9[) = S+ido,9[)\S+{do,9i) 



then 



= 1 - e--^l^+\-S+(rfo,ei)| ^ ^ 



(12) 



It is the complement of the probability that no node is in 
the area A+\S+{da,9i). 

From (|8j, (|9|, ([TTJ and ([T2|, we can obtain 

P{ U T((ro,^),Ti,n2)} 

+ P{ y T((ro,7r),T(,n2) I r'((ro,7r),Ti,n3)} 

{n2.r>3 }£«■!.(, n 



loiro) 



(13) 



IEEE/ACM TRANSACTIONS ON NETWORKING 



7 



Similarly, we can further consider the case that neither of ti and 
and t[ can contribute to a triangle which bounds a triangular 
hole containing the origin, but other nodes with even higher Pu{X,j) = 2tt\^ 
polar angle can contribute to such a triangle. In this way, we 



can get the same result as ( 13 i. 
Therefore, it can be derived that 

P{ y r((ro,^),ni,n2)} 

{ni,«2}C$' 



< 



D(ro) 



(14) 



From (j2]i, (jSj), (j6]l, (|7]) and ( 14 1, the upper bound shown in 
Q can be derived. 

Here we need to compute the areas of ^^(ro,^!), 
5^(ro,ri,6'i) and S'^(ro, rg, 6'i). In fact, the areas 
|5+(ro,0i)|, |5-(ro,ri,0i)| and |5-(ro, tq, have 
very similar expressions. For example, the area |5+(ro,0i)| 
can be expressed as 

|5+(ro,0i)| = / dO rdr 

= 1 [ \l{ro,e,)d9-'fi9,^ao) 
When 6*1 < 7r/2 - arccos(i?c/(2ro)) 

|5+(ro,0i)| =/(0i)-/(ao) 

where 



, 7'o sm f 



m\ ==-0 sin 61 cos 61 + arcsin 

^ ' 2 2 2 Rr 



SVa9\ Wi — Tn SIT? 



2 V 2 
When 7r/2 - arccos(i?c/(2ro)) < 6*1 < ai 

|^+(ro,0i)| =2|5+(ro,7r/2-arccos(i?,/(2ro))| 

- |^+(ro,7r - 2arccos(i?,/(2ro)) - 0i| 

Similarly, |S'^(ro, ri, 6*1)! and |S'^(ro, ro, 6'i)| can be ob- 
tained. ■ 



D. Case 7 > 2 

Theorem 3. When 7 > 2, pi{X,j) < p(A,7) < p„(A,7), 
w/iere 

Pi{X,j) = 2ttX^ / rodro / d0i / 

^ Ji?,., Jq Jra 

Ra/V^ pai /--Rl(r-0:ei) ^^^^ 

rodro / d^i 



-A|S+(ro,ei)| 



(l_e-^l^"(''«^'^i>«i)l)ridri} 



rodro / rf6'i / 

Rs Jo Jro 

-\\S+{ro, 01)1(^1 _ e-MS-iro.ro,e,)\y^^^^ 

Ra/VS i-ai /■-Ri(ro,ei) 

„-Airr„ 



rodro / d6li 

-Rc/2 -^Qo -Jro 

+ P{ (J r(no,ni,n2) I T'(To,n3,n4)} 

{no,--- ,"4}C*\{to($)} 

(16) 



and 



Ri{ro, 0i) = mm^J R^ - sin 9i — tq cos^i, 



i?2 - sin^ 01 + ro cos 6*1) 




Fig. 6. Illustrations of areas in case 7 > 2. (a) the areas j4+, S"*" and 5 
(b) the areas S'+ and S'" 

In this case, we can use the same method as in Section 
IV.C to get the lower and upper bounds, shown in ( fTSj ) and 
(16 1 respectively. But we need to consider two situations Rg < 
dn < Rc/2 and Rc/2 < dg < Rc/VS. In the first situation, 
do < Rc/2 means that the ball B{0, do) is included in the ball 



g 
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B{to,Rc). The illustrations for the areas 5+, , 5'+ and 
5'^ are shown in Figure|6|a) and|6|b) respectively. In addition, 
the lower limit of integration for 9i is and the upper limit 
is TT. The second situation is the same as that in Section IV.C. 

V. Distributed coverage hole detection 

ALGORITHM 

Our algorithm includes five components: weight computa- 
tion, vertex and edge deletion, boundary edge detection, coarse 
boundary cycles discovery and boundary cycles minimization, 
as shown in Figure |7] An example is used to illustrate the 
procedures of this algorithm in Figure |8] For a WSN with 
some coverage holes, such as the one shown in FigurejHJa), our 
algorithm aims to discover minimum boundary cycles of all 
holes. In weight computation component, each node computes 
its weight independently. The definition of weight of a node 
will be presented in the next part. After obtaining the weight, 
each node continues to determine whether it can be deleted 
or not according to some rules defined hereafter. Figure |8jb) 
shows the result of vertex deletion. Furthermore, some special 
node will decide whether some edge can be deleted or not. 
Figure [Hj^c) shows the process of such special edge deletion. 
After the second component, many boundary edges can be 
found, as the red line shown in Figure [8|d). But it is possible 
that some other boundary edges have not been found. Then 
in the third component, all or nearly all boundary edges will 
be found after deleting some edges, see Figure |8je-j). Then 
coarse boundary cycles can be easily discovered, as shown in 
Figure [8|k). It is possible that the found boundary cycles are 
not minimum. In this case, coarse boundary cycles will be 
minimized in the final component as shown in Figure [Hj^l). 



Start 



Weight 
computation 

r 

Vertex and edge 
deletion 

T 

Boundary edge 
detection 



Coarse boundary 
cycles discovery 

I 

Boundary cycles 
minimization 
I 
End 



Fig. 7. Flow chart of the algorithm 



A. Definitions 

Before presenting the details of our algorithm, we first 
give some definitions that will be used in the process of the 
algorithm. 

We say that a i-simplex [vio, vn, Vii]is pait of a j-simplex 



[vjo,Vji, ...,Vjj] if [vio.Vii, ...,vu] C [vjo,Vji, 



"33 



]. So the 



vertex vq and vi is part of the edge vaVi. The edge vovi is 
part of the triangle Vf)ViV2- In addition, we use E{v) to denote 
all the edges that the node v is part of and T{y) to denote all 
the triangles that the node v is part of. 



Definition 4 (Index of a triangle). The index of a triangle A 
is the highest dimension of the simplex that the triangle is part 
of, denoted by 

Definition 5 (Weight of a node). The weight of a fence node 
is defined to be 0. For any internal node v, if there exists one 
edge in E{v) which is not part of any triangle, the weight 
of node v is set to be 0; if not, the weight is the minimum 
index of all the triangles in T{v), that is Wy — min^^Tiv) I A 

The weight of an internal node is an indicator of the density 
of surrounding nodes. If the weight of an internal node is 0, 
the node must be on the boundary of a coverage hole. The 
larger the weight is, the higher the probability that the node 
is not on the boundary of a coverage hole. 

We also use the definition of simple-connectedness graph 



as in |27|. Let G be a simple graph with vertex set V{G) 
and edge set E{G). A cycle C is a sub-graph of G if it is 
connected and each vertex in C has degree two. The length 
of a cycle C is the number of its edges, denoted by |i?(C)|. 
The cycle space G{G) of a graph G contains all the cycles 
in G. The addition of two cycles Ci and C2 is defined as 
Ci © C2 = {E{Ci) U E{C2)) \ E{Ci) n E{G2). The triangle 
cycle subspace Ct{G) of G is the set of all 3-length cycles 
in C(G). 

Definition 6 (Simple-Connectedness Graph). A connected 
graph G is of simple connectedness if its cycle space G (G) 
is empty, or for any cycle G in C{G), there exists a set of 
3-length cycles Tq C Ct{G) such that G = X^tgTo 

Let X be a vertex (or edge) set in a graph G, we use G[X] to 
denote the vertex-induced (or edge-induced) sub-graph by X. 
The neighbour set of a vertex w in G is denoted by Na{v). The 
neighbouring graph rQ{v) of vertex v is denoted as G[Nq{v)]. 
The neighbouring graph Tc{e) of an edge e — {u, v) is defined 
as G[Nc{u) n Ng{v) U — e. The neighbour set of k- 

simplex [vo,vi, ...,Vk] is defined as PliLo ^G('f^i)- 

Definition 7 (Deletion of fc-simplex in Rips complex TZ{V)). 
A k-simplex [vq,Vi, ■ ■ ■ , w^] is deleted in a Rips complex TZ{V) 
means that the simplex and all simplices which the simplex is 
part of are deleted from TZn^iV). 

Based on definitions above, we can give the definition of 
HP (Homology Preserving) transformation. 

Definition 8 (HP Transformation). A HP transformation is a 
sequential combination of vertex (or edge) deletion as follows: 
a vertex (or edge) x of G is deletable if neighbouring graph 
Tq{x) (1) has two or more vertices; (2) is connected and (3) 
is a simple-connectedness graph. 

Tlieorem 4. HP transformations do not change the number 
of coverage holes in Rips complex of a WSN. 

Proof: In order to prove HP transformations do not 
change the number of coverage holes in Rips complex of 
a WSN, we only need to prove that in the process of any 
HP transformation, there is no new coverage holes created 
and no two coverage holes merged. If a new coverage hole 
is created when a vertex v (or edge e) is deleted, then the 
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boundary cycle of the new coverage hole must be a cycle 
in Tg{v) (or rG(e)), which means Tg{v) is not a simple- 
connectedness graph. It is contrary to the third condition in 
HP transformation, so there is no new coverage hole created. 
Furthermore, if two coverage holes are merged when a vertex 
V (or edge e) is deleted, then the neighbour graph Tq{v) 
(rG(e)) must not be connected, which is contrary to the second 
condition in HP transformation. So no two coverage holes are 



merged. Consequently, the number of coverage holes will not 
be changed in the process of any HP transformation. ■ 

B. Weight computation 

In this component, each node computes its weight. For any 
fence node, its weight is 0. For any internal node, theoretically 
it needs to construct all the simplice it is part of. As we 
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consider WSNs in a planar target field, each internal node only 
needs to construct all its 1 -simplex and 2-simplex and their 
neighbours. This can also reduce the computation complexity. 
In order to do this, the node needs to obtain all its 1- and 
2-hop neighbours information. This can be easily achieved 
by two broadcasts of hello message. In the first one, each 
node broadcasts its id. When it gets all the ids of its 1-hop 
neighbours, each node continues to broadcast a hello message 
containing the ids of its 1-hop neighbours. After receiving the 
neighbour list of its neighbours, the node can obtain its E{v), 
the set of edges (1 -simplex) and T{v), the set of triangles (2- 
simplex). It can also get the neighbours of each simplex. For 
any e e E{v), let n{e) denote the neighbour set of e. For any 
t e T{v), let n{t) denote the neighbour set of t. Then the 
weight of node v can be computed as in Algorithm [T] 

Algorithm 1 Weight computation (for internal node v) 
Begin 

if 3e G E{v),n{e) is empty then 

Wy = 

else if 3t e T{v),n{t) is empty then 

Wy = 2 

else 

Wy = 3 

end if 

END 



C. Vertex and edge deletion 

In this component, we conduct maximal vertex deletion 
without changing the number of coverage holes in the original 
WSN and also delete some special edges if such edges exist. 

1) Vertex deletion: From the definition of weight, we can 
see that the higher the weight is, the higher is the probability 
that the sensing range of the node is fully covered by its 
neighbours, consequently the probability that the node does 
not he on the boundary is higher. Meanwhile, if the deletion 
of a vertex may create a new coverage hole, it must not be 
deleted no matter how high the weight is. So we have such a 
rule for vertex deletion. If the weight of a vertex is smaller than 
3, it should never be deleted. Otherwise, the vertex continues 
to check whether it is deletable or not according to HP 
transformation. After the verification, the vertex broadcasts a 
message indicating that it can be deleted or not. After receiving 
the status of all its neighbours, each deletable vertex continues 
to check whether it should be deleted. It can be found that the 
weight of any deletable vertex must be 3. We assume that the 
vertex with lower ID has the priority to be deleted first. So 
each deletable vertex just needs to check whether its ID is 
the lowest among all its deletable neighbours. If so, it should 
be deleted. Otherwise, it should not be deleted. Algorithm |2] 
gives the detailed process for vertex deletion. According to 
the rule, it can be seen that two neighbouring vertices will 
not be deleted simultaneously, so each vertex can make the 
decision independently. When a vertex is deleted, it broadcasts 
a message to its neighbours. All its neighbours will modify 
their simplices accordingly and compute their weights again. 



The procedure of vertex deletion terminates until no vertex 
can be deleted in the Rips complex. Figure [8]^b) gives the 
final result after vertex deletion. 



Algorithm 2 Vertex deletion (for internal node v) 
Begin 

if lUt, < 3 then 

node V can not be deleted 
else if node v is not deletable according to HP transforma- 
tion then 

node V can not be deleted 
else if the ID of node v is the smallest among all its deletable 
neighbours then 

node V is deleted 
end if 
END 



2} Edge deletion: After vertex deletion, we find one inter- 
esting thing. It is highly possible that an edge having only one 
neighbour lies on the boundary of a coverage hole, such as the 
red edges shown in Figure |8jc). It is also possible that there 
exist some special such edges not lying on the boundary, such 
as the blue edges shown in Figure [8]^c). We try to delete such 
special edges. Since the edge has only one neighbour, deleting 
the edge will not change the number of coverage holes. We 
design a rule for edge deletion. If a vertex has only one special 
edge and deletion of the special edge can make the vertex have 
a Hamilton cycle in its left neighbouring graph, then such a 
special edge can be deleted. According to this rule, nearly all 
edges which have only one neighbour and are not lying on the 
boundary of any holes can be deleted. But it is also possible 
that some such edges lying on the boundary are also deleted. 
This is not a big issue, because deletion of such edges will 
not change the number of coverage holes and just enlarge the 
current coverage holes. It can be solved in the boundary cycles 
minimization component. In addition, after edge deletion, it 
is possible that some vertices can be deleted again, such as 
the vertex denoted by red square in Figure |8jc). If such a 
case happens, we can continue to do vertex deletion until no 
more vertex or edge can be deleted. Figure [SJd) shows the 
result after edge deletion. Here we need to talk a little more 
about Hamilton cycle. In general, it is a NP-complete problem 
to find a Hamilton cycle in a graph. However, considering 
the left nodes after vertex deletion usually have very small 
number of neighbouring nodes, it is not complicated to check 
the existence of a Hamilton cycle in their neighbouring graph 
even if we use exhaustive search method. 

D. Boundary edge detection 

After deleting some vertices and edges, we continue to look 
for the edges which have only one neighbour It can be seen 
that nearly all such edges lie on the boundary of holes. We 
call them boundary edges. It can also be found that some 
edges lying on the boundary have not been found. In this 
component, we try to find such edges as many as possible. 
In all cases, such edges have two or more neighbours. If we 
consider the nodes having one or more boundary edges as 
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boundary nodes and other nodes as non-boundary nodes, then 
we try to delete some edges connecting non-boundary nodes 
and boundary nodes according to HP transformation, such as 
the green edges shown in Figure |8je). After that, it is possible 
that some new boundary edges are recognized as shown in 
Figure [Hj^f). And it is also possible that some edges lying 
on the boundary have still not been discovered. Such case 
usually occurs when some boundary nodes are neighbours and 
edges connecting them have more than one neighbours. In this 
case, we randomly delete some of such edges according to HP 
transformation, such as the green edges in Figure |8jg). After 
that, some new boundary edges can be recognized, as shown 
in Figure |8|h). 

But it is possible that the new found boundary edges are 
not correct, as the magenta edge in coverage hole 5 in Figure 
|8jh), or can not construct a cycle with other boundary edges, 
as the magenta edges in coverage hole 1, 2 and 3 in Figure 
|8jh). As for the case in coverage hole 5, we define such a 
rule to delete a boundary edge. For any internal node v, if it 
has two boundary neighbour nodes u and w, the edge vu and 
vw are not boundary edge and wu is a boundary edge, if the 
deletion of the edge wu can make at least one of the two edges 
vu and vw be boundary edge, then wu can be deleted. As for 
the case in coverage holes 1, 2 and 3, it is often due to the 
fact that two boundary edges cross with each other. It is found 
from numerous simulations that there are mainly three cases 
of crossing boundary edges, as shown in Figure |9j red lines 
denote boundary edges and black ones denote non-boundary 
edges. Similarly, we define some rules to delete such edges. 
We still take node v as an example, in the top part of Figure 
[9]]a), the two red lines connecting v and its neighbours are 
deleted, as shown in the bottom part of Figure |9|a). For the 
cases in Figure |9|b) and (c), the red line connecting v and its 
neighbour is deleted, as shown in the bottom part of Figure 
[9]^b) and (c). According to such rules, some boundary edges 
can be deleted, such as the black bold edges in Figure |8ji). 
From Figure [8]^i), it can be found that certain boundary edge 
is deleted incorrectly. It is not a big issue as explained in 
last part. Deleting such edges will not change the number of 
coverage holes and just enlarge the current holes, it will be 
solved in the process of minimizing boundary cycles. After 
deletion of such edges, new boundary edges can be found as 
shown in Figure IHl 



\:><7 ^xl 
'4 "I ' 



(a) (b) (c) 

Fig. 9. Illustration of crossing boundary edges 

In general, after the process above, nearly all boundary 
edges can be found. But there exists one special case as in 




Figure 10 In this case, some edges lying on the boundary can 
not be detected. It can be solved in the next component. 



Fig. 10. Special case when some boundary edges can not be detected 

E. Coarse boundary cycles discovery 

After boundary edges are detected in the former component, 
it is easy to discover the coarse boundary cycles. We just need 
to randomly choose one node which has two boundary edges 
in any boundary cycle. The node initiates the process to find 
the coarse boundary cycle by sending a message along one of 
the boundary edges. Then the boundary neighbour continues 
sending the message along its boundary edges. When the 
initiating node receives the message coming back along the 
other boundary edge, it discovers one coarse boundary cycle. 
Similarly, all coarse boundary cycles can be found, as shown 
in Figure |8jk). 

As for the special case shown in Figure [TO] when the node 
V receives a message from its boundary neighbour node u, it 
broadcasts the message to all its neighbours except u. If its 
neighbour node is a boundary node, then the message can be 
sent along the boundary edges. If its neighbour node is not a 
boundary node but it has boundary neighbour nodes, then it 
can send the message to its boundary neighbour nodes. Else, it 
will not transmit the message again. In this way, the message 
goes along boundary edges most of the time and can return 
to the original node initiating the message. 

F. Boundary cycles minimization 

It is possible that some primary boundary cycles found are 
not minimum, so we need to minimize such cycles. This can 
be achieved by checking whether there exists shorter path 
between any two nodes in the cycle. Since each node has 
its 1- and 2-hop neighbours information, it can locally check 
the existence of shorter path in the cycle in general cases. 
If there exists, we shorten the cycle and continue to do the 
same verification until no such case exists. After that, it is 
still possible some cycle has not been minimized, such as the 
coverage hole 2 in Figure [HJk). So we use the following 2- 
hop shrinking process to make the cycle as shorter as possible. 
For any four adjacent nodes in the cycle, say a, 6, c, d, if there 
exists one node x which is one common neighbour of nodes 
a, 6, c, d, then the cycle can be shortened by using x to replace 
nodes b and c. 

In this way, we can nearly obtain most minimum cycles 
surrounding coverage holes. It is also possible that in some 
cases, we can not get the minimum cycles since each node 
only has its I- and 2-hop neighbours information. Even so, the 
boundary cycles discovered in the algorithm can still provide 
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valuable information about coverage holes. Theoretically, if 
each node has more information about its fc-hop neighbours, 
aU cycles can be minimized. 

VI. Simulations and performance evaluation 

Performance evaluation of the theoretical bounds obtained 
in Section IV and the algorithm proposed in Section V is 
presented in this section. 



extent, it is unacceptable to use connectivity information only 
to discover coverage holes. 



Finally, it can be found in Figure 1 1 a) that the probability 
obtained by simulation is very well consistent with the lower 
bound. The maximum difference between them is about 0.5%. 
Figure [TTJb) shows that probability obtained by simulation is 
also consistent with the upper bound. The maximum difference 
between them is about 3%. 



A. Simulation settings 

We have different simulation settings to evaluate the bounds 
and the algorithm. 

For bounds evaluation, a disk centered at the origin with 
radius Rc is considered in the simulations. The probability 
that the origin is inside a triangular hole is computed. Sensors 
are randomly distributed in the disk according to a Poisson 
point process with intensity A. The sensing radius Rg of each 
node is set to be 10 meters and 7 is chosen from 2 to 3 
with interval of 0.2. So the communication radius Rc ranges 
from 20 to 30 meters with interval of 2 meters. A is selected 
from 0.001 to 0.020 with interval of 0.001. For each 7, 10^ 
simulations are run under each A to check whether the origin 
belongs to a triangular hole. 

For performance evaluation of the algorithm, we choose a 
100 X 100 m^ square area as the target field. The sensing 
radius Rg of each node is 10 meters. The communication 
radius Rc is set to be 20 meters and so 7 = 2. There are 
fence sensors locating along the edges of the square with 20 
meters distance between neighbours. Other internal sensors are 
randomly distributed in the area according to a Poisson point 
process with intensity A. 

B. Proportion of the area of triangular holes 

The probability p(A, 7) obtained by simulation is presented 
with the lower and upper bounds in Figure [TTJa) and [TT|b) 
respectively. It can be seen that for any value of 7, p(A, 7) 
has a maximum at a threshold value Ac of the intensity. 

As a matter of fact, for A < Ac, the number of nodes is 
small. Consequently the probability that the origin belongs to 
a triangular hole is relatively small too. With the increase of A, 
the connectivity between nodes becomes stronger As a result, 
the probability that the origin belongs to a triangular hole 
increases. However, when the intensity reaches the threshold 
value, the origin is covered with maximum probability. p(A, 7) 
decreases for A > Ac. The simulations also show that Ac 
decreases with the increase of 7. 

On the other hand, it can be seen from Figure [TTfa) and 
1 1 b) that for a fixed intensity A, p(A, 7) increases with the 
increases of 7. That is because when Rs is fixed, the larger 
Rc is, the higher is the probability of each triangle containing 
a coverage hole. 

Furthermore, the maximum probability increases quickly 
with 7 ranging from 2.0 to 3.0. It is shown that when 7 = 2, 
the maximum probability from simulation is about 0.03% and 
thus it is acceptable to use Rips complex based algorithms to 
discover coverage holes. While the ratio 7 is high to a certain 
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Fig. 1 1 . Proportion of the area of triangular lioles (a) simulation results and 
lower bounds ; (b) simulation results and upper bounds 

C. Performance of the algorithm 

1 ) Complexity: The computation complexity of each step in 
the algorithm is shown in Table |l] In the weight computation 
component, each node only needs to check all its 2-simplex, so 
the computation complexity is 0{n'^), where n is the number 
of its 1-hop neighbours. 

In vertex deletion part, each node needs to check whether 
it is deletable or not according to HP transformation. This 
can be done by checking all its cycles in its neighbouring 
graph. It can build a spanning tree in its neighbouring graph 
and check all fundamental cycles in the spanning tree. There 
are E — n + 1 fundamental cycles, where E is the number of 
edges in its neighbouring graph, so the worst case computation 
complexity is O(n^). Since the node needs to recompute its 
weight and recheck whether it is deletable when any one of 
its neighbour is deleted, so the total worst case computation 
complexity is 0{n^). As for the edge deletion, only limited 
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TABLE I 

COMPLEXITY OF EACH STEP IN THE ALGORITHM 



Step 


Complexity 


Weight computation 


O(n^) 


Veitex and edge deletion 


O(n^) 


Boundary edge detection 


O(n^) 


Coarse boundary cycles discovery 


0(1) 


Boundary cycles minimization 


0(1) 



nodes will do this process and the complexity is not high as 
explained in Section V.C. 

In the boundary edge detection component, the non- 
boundary nodes need to check whether the edge connecting 
itself with its boundary neighbours can be deleted or not 
according to HP transformation. So the worst case computation 
complexity is 0{n^), the actual complexity is much less than 
that since for one edge, there are usually very few nodes 
in its neighbouring graph. The boundary nodes also need to 
check whether the edges connecting itself and its boundary 
neighbours can be deleted or not, so the complexity is also 
0{n^). In addition, the boundary nodes need to check whether 
there exist special cases as illustrated in Section V.D. The 
worst case computation complexity for such process is 0{n^). 
So the complexity of this step is 0{n^). 

As for the final two components, each node only needs to 
broadcast some messages and do some local computations, 
the complexity is 0(1). So the total worst case computation 
complexity for our algorithm is 0{n^). 

2) Comparison with other algorithms: In order to evaluate 
the performance of our proposed homology based algorithm 
(denoted as HBA), we compare it with the location based 
algorithm (denoted as LBA) proposed in |19| . Since loca- 
tion based algorithm can discover both triangular and non- 
triangular coverage holes, but our algorithm can only detect 
non-triangular coverage holes, we do not consider those tri- 
angular coverage holes in the comparison. Since it is possible 
that there exist shorter paths in boundary cycles found by 
LBA, we first shrink them using 1-hop neighbour information 
of boundary nodes. After that, we compare those boundary 
cycles with what our algorithm finds. Since for some coverage 
holes, the minimum boundary cycles may not be unique, two 
boundary cycles are considered to surround the same coverage 
hole if one cycle can be converted to another one by using 
only 1-hop neighbours information. We emphasize that only 
1-hop neighbours information can be used in the comparison in 
order to evaluate the accuracy of boundary cycles found by our 
algorithm. For example, if one cycle ci found by our algorithm 
can not be converted to another cycle C2 found by LBA using 
only 1-hop neighbours information but can be converted by 
using 2-hop neighbours information, we consider the cycle ci 
is not accurate and the corresponding coverage hole is not 
found. 

Based on the method presented above, we set A to be 
0.008, 0.010 and 0.012 to represent sparse, moderate and dense 
WSNs respectively. For each intensity, 1000 simulations are 
performed. Simulation results show that when A is 0.008, there 
are nine times among the 1000 times when our algorithm can 
not find all non-triangular coverage holes. In each of the nine 



times, only one coverage hole is missed. There are 7363 non- 
triangular holes in total and 7354 ones found by our algorithm. 
When A is 0.010 and 0.012, only one time among the 1000 
times when our algorithm can not find all coverage holes. 
And in that time, only one coverage hole is missed. When 
A is 0.010, there are 6114 non-triangular holes in total and 
6113 ones found by our algorithm. When A is 0.012, there 
are 4613 non-triangular holes in total, of which 4612 ones are 



found. The results are shown in Figure 12 All these results 
show that our algorithm can find about 99% coverage holes 
in about 99% cases. 
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Fig. 12. Performance evaluation of the algorithm 



VII. Conclusions 

In this paper, we adopt two types of simplicial complex 
called Cech complex and Rips complex to capture coverage 
holes of a WSN. The relationship between Cech complex 
and Rips complex in terms of coverage hole is first analysed 
under different ratios between communication radius and 
sensing radius of a sensor Based on that, we define two 
types of coverage holes: triangular and non-triangular hole. 
For triangular holes, both the lower and upper bounds on 
the proportion of the area of triangular holes in a WSN 
are derived. Such proportion is related to the ratio between 
communication radius and sensing radius of each sensor When 
the ratio is no larger than ^/3, there is no triangular hole. When 
the ratio is between \/3 and 2, both the theoretical analysis 
and simulation results show that the proportion is lower than 
0.06% under any intensity. It means that the triangular holes 
can nearly be neglected. When the ratio is larger than 2, the 
proportion of the area of triangular holes increases with 7. It 
becomes unacceptable for 7 larger than a threshold. In that 
case triangular holes can not be neglected any more. For non- 
triangular holes, a homology-based algorithm is proposed to 
detect them. Simulation results show that the algorithm can 
detect 99% such holes. 
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